Phase Ordering and Onset of Collective Behavior in Chaotic Coupled Map Lattices 
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The phase ordering properties of lattices of band-chaotic maps coupled diffusively with some 
coupling strength g are studied in order to determine the limit value <? e beyond which multistability 
disappears and non-trivial collective behavior is observed. The persistence of equivalent discrete 
spin variables and the characteristic length of the patterns observed scale algebraically with time 
during phase ordering. The associated exponents vary continuously with g but remain proportional 
to each other, with a ratio close to that of the time-dependent Ginzburg-Landau equation. The 
corresponding individual values seem to be recovered in the space-continuous limit. 



One of the most remarkable features distinguishing 
extensively-chaotic dynamical systems from most mod- 
els studied in out-of-equilibrium statistical physics is 
that they generically exhibit non-trivial collective be- 
havior (NTCB), i.e. long-range order emerging out of 
local chaos, accompanied by the temporal evolution of 
spatially-averaged quantities Jl|-§|. In particular, NTCB 
is easily observed on simple models of reaction-diffusion 
systems such as coupled map lattices (CMLs) in which 
(chaotic) nonlinear maps S of real variables X are cou- 
pled diffusively with some coupling strength g || . 

NTCB is often claimed to be a macroscopic attractor, 
well-defined in the infinite-size limit and reached for al- 
most every initial condition, provided the local coupling 
between sites is "large enough" . On the other hand, for 
small g values, such as those corresponding to the so- 
called "anti-integrable" limit which tries to extend zero- 
coupling behavior to small, but finite coupling strengths, 
CMLs often exhibit multistability ||. This is in par- 
ticular the case if the local map shows banded chaos, 
because the interfaces separating clusters of sites in the 
different bands can be pinned. This multistability is "ex- 
tensive" : the number of (chaotic) attractors may then 
be argued to grow exponentially with the system size, in 
opposition to NTCB for which this number is small and 
size- independent . 

In this Letter, we define and measure the limit cou- 
pling strength g c separating the strong-coupling regime 
in which NTCB is observed from the weak-coupling, ex- 
tensive multistability region. Using the discrete "spin" 
variables which can be defined whenever the one-body 
probability distribution functions (pdfs) of local (contin- 
uous) variables have disjoint supports, we study numer- 
ically the phase ordering process following uncorrelated 
initial conditions in cases where the spin variables take 
only two values. We find that the persistence probability 
p{t) (i.e. the proportion of spins which have not changed 
sign up to time t) saturates in finite time to strictly pos- 
itive values in the weak coupling regime whereas it de- 



cays algebraically to zero when g > g c . The associated 
persistence exponent 9 varies continuously with parame- 
ters, at odds with traditional models ||. Moreover, data 
obtained on various two-dimensional CMLs is best ac- 
counted for by a relation of the form 9 ~ (g — g e ) w , 
which we use to estimate g e . We show further that this 
behavior is mostly due to the non-trivial scaling of the 
characteristic length Lit) ~ t^ during the phase ordering 
process. Indeed, <j> ^ \, the expected value for a scalar, 
non-conserved order parameter H, and is found to be 
proportional to 9, with the exponent ratio <p/6 approx- 
imately taking the value known for the time-dependent 
Ginzburg-Landau equation (TDGLE). We also provide 
evidence that, in the continuous-space limit, "normal" 
phase ordering behavior is recovered. Finally, we discuss 
the hierarchy of limit coupling values <?" which can be 
defined when the local map is unimodal and shows 2™- 
band chaos, using recent results on renormalisation group 
(RG) ideas applied to CMLs @. 

Consider a d-dimensional hypercubic lattices C of cou- 
pled identical maps Sn acting on real variables iXp)p e c : 

A< +1 = (1 - 2dg)S tl {Xi) + g £ S^Xi +s ) , (1) 

where V is the set of 2d nearest neighbors e of site 0. 
We first present results obtained for the piecewise linear, 
odd, local map defined by: 

( iiX if X € [-1/3,1/3] 

S„{X) = { 2^/3-AiX if Xe [1/3,1] (2) 
[ -2/z/3 - fiX if X e [-1, -1/3] 

which leaves the / = [—1, 1] interval invariant. (For 
(i = 3, this is the chaotic map introduced by Miller and 
Huse M.) For fi e [—2,-1], displays banded chaos, 
while for opposite fj, values, these bands become invariant 
subintervals of /. At /i = 1.9 in particular, possesses 
two symmetric such intervals I — [±/i(2 — //)/3,±/i/3], 
separated by a finite gap. For any value of g, the sup- 
port of the pdf of X for the CML defined by (§f|) can be 
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separated into two components thanks to the symmetry 
of the map. This allows the unambiguous definition of 
spin variables oy = sign (.XV). The deterministic nature 
of the system and the form of the coupling strictly forbids 
the nucleation of opposite-phase droplets in clusters: the 
analog spin system is at zero temperature. 




FIG. 1. Snapshots of the d = 2 CML with local map (g). 
Lattice of 128 2 sites, grey scale from X = — 1 (white) to X = 1 
(black), uncorrelated initial conditions. (a,b): transient lead- 
ing to complete ordering at g = 0.2 > g e , t = 100 and 1000; 
(c,d): blocked state at g = 0.15 < g e , t = 1000 and 2000. 
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FIG. 2. Phase ordering in the d = 2 CML with local map 
(^) at fj, = 1.9. (a): algebraic decay of the persistence proba- 
bility p(t) for, from top to bottom: g — 0.16 (< g c : saturates 
to a finite level), 0.17, 0.172, 0.174, 0.176, 0.18, 0.185, and 
0.195. (b): variation of persistence exponent with g. Solid 
line: fitting Ansatz ~ (g — g e ) w with g c = 0.169(1) and 
w = 0.20(3) Insert: log(0) vs. log(g - g e ). 

For large g values, complete phase ordering occurs 
(Fig. ^a,b), and the system eventually reaches a regime in 
which all sites are situated in one of the two intervals 1^ . 
For small g, initial conditions with sites in both intervals 
1^ lead to spatially-blocked configurations where inter- 



faces between clusters of each phase are strictly pinned, 
while chaos is present within clusters (Fig. |l|c,d). (!(]] 

To study the phase ordering process efficiently, uncor- 
related initial conditions were generated as follows: ex- 
actly one half of the sites of a d = 2 lattice were chosen at 
random and assigned positive X values drawn according 
to the invariant distribution of on I + , while the other 
sites were similarly assigned negative values. Large lat- 
tices with periodic boundary conditions were used, and 
the persistence pit) was measured. Fig. |]a shows the 
results of single runs for various values of g. For small 
g, p(t) saturates at large times to strictly positive values, 
while it decays algebraically, for large g, on square lattices 
of linear size 2048 sites. The associated persistence expo- 
nent 9 varies continuously with g, and its ^-dependence 
is nicely accounted for by a functional form 9 ~ (g — g c ) w 
with g c ~ 0.169(1) and w ~ 0.20(3) (Fig. |b). We have, 
at this point, no theoretical justification of this fitting 
Ansatz. At any rate, it provides an operational definition 
of g c yielding estimates consistent with those obtained 
using other, less accurate, methods [ O . 



(a) (b) 




7 2 3 4 0.2 0.3 0.4 



log* 4> 

FIG. 3. Same runs as in Fig. pi (a) log(L) vs log(£) (bot- 
tom curve: g = 0.16, top curve g = 0.195); (b) 0(g) vs <t>{g) 
for g values between 0.17 and 0.20; the solid line is the linear 
fit ~ 0.3960 - 0.002. 

The origin of this unusual behavior of the persistence 
exponent is largely explained by the evolution of the spa- 
tial structures formed during phase ordering. Usually, 
one expects the coarsening to be described by the alge- 
braic growth of a single characteristic length L(t) ~ If 
with 0=1/2 for a non-conserved, scalar order parameter 
[Q. In the CML studied above, the two-point correlation 
function C(x, t) = (a~ +s ap) was measured during phase 
ordering Length L(t) was then evaluated to be the 
width at mid-height (C(L(t),t) = 1/2), determined by 
interpolation. This procedure was then validated by a 
collapse of all C(\\x\\/L(t),t) curves. Surprisingly, while 
the scaling behavior of Lit) is observed, exponent </> de- 
parts from the expected 1/2 value and varies continu- 
ously with g (Fig. ^). Again, we find a law of the form 
<j) ^ (g—g a ) w to be an acceptable Ansatz of our numerical 
results. The estimated values of g e and w are consistent, 
within numerical accuracy, with those found when fitting 
9(g). This is corroborated by studying directly p(t) vs 
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L(t) (not shown), or by plotting 9 vs cf> which confirms 

that the two exponents are proportional to each other 

(Fig. ||d). Remarkably, the ratio 9/(f> is found to have, 

within our numerical accuracy, the d = 2 TDGLE value: 

9/4> ~ 0.40(2) ~ 29 GL ^ 0.40 Q. (We cannot, however, 

completely exclude the values corresponding to the Ising 

model, or the diffusion equation, since #i S ing — 0.22 Jli) , 

and fciff. ^ 0.19 O.) 

(a) ^ (b) 
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FIG. 4. d = 2 lattice of coupled logistic maps for g = 0.2. 
(a) bifurcation diagram ((Xp)?, large dots) superimposed on 
that of single logistic map (small dots), (b): (er£)p vs p plot- 
ted at 20 different timesteps between t — 500 and t = 1000 
for < /tt = 1.5 < fiz. This "magnetization" remains con- 
stant only for p = p* ~ 0.3465, whereas it reflects the overall 
period-2 dynamics for all other p values, (c) 6(g) at /i = 1.5 
and p = p*(g) with log-log fit (g e ^ 0.101, tu ~ 0.06). (d) 
0(g) vs 0(g); the solid line is the linear fit 9 ~ 0.400 - 0.002. 

The same analysis was also performed on CMLs with 
a non-symmetric, unimodal, local map of the form: 

= 1 - p\X\ 1+£ with fx e [0, 2] , (3) 

in particular for e = (tent map) and e = 1 (logis- 
tic map) . For p S [/ioo , 2] , this map shows 2"-band 
chaos and exhibits an inverse cascade of band-merging 
points /Et„ when p — * /ioo. In the strong-coupling limit, 
the corresponding CMLs exhibit, depending on d, peri- 
odic or quasiperiodic NTCB with a period equal to, or 
a multiple of, that of the band-chaos of the local map 
HD- For d = 2 and 3, in particular, simple period-2 n 
NTCB occurs, with an infinite cascade of phase tran- 
sition points p c n distinct from the band-merging points 
(Fig. ||a). When period-2 NTCB occurs in the two-band 
chaotic region of the map (pi S [p,2,pi] ~ [1.43,1.54]), 
two-state spin variables at g { — 1, 1} can be defined, but 
the asymmetry of the two bands hinders the generation 
of "effectively" uncorrelated initial conditions. Indeed, 
an equal proportion of sites in each band quickly leads to 
complete phase ordering and saturation of p(t), even in 
the strong-coupling regime. This happens because these 
initial conditions create, after a few timesteps, config- 
urations with a fairly large unbalance between the two 
phases. Tuning the initial proportion p of sites in, say, 
the band containing X = 0, one can minimize such ef- 
fects. We determined the optimal proportion p* defined 
as the value for which the magnetization (at) remains 



constant (Fig. ||b). Clean scaling behavior of L(t) and 
p(i) is then observed with reasonable system sizes, as 
with the symmetric local map (|J). Varying the coupling 
strength <?, exponents <p and 9 show the same behavior 
as above, decreasing continuously to zero at g c . Fig. [|c 
shows the case of coupled logistic maps, for which the 
Ansatz 9,<fr ~ (.9 — g c ) w is, again, valid, although not as 
good as in the case of local map (||) . Note that the esti- 
mated value w ~ 0.06(2) is different from that measured 
for the CML with local map (§), but 6/<j> ~ 0.48(4) is 
still rather close to the TDGLE value (Fig. |d). 

We now deal with the onset of more complex NTCB 
such as the period-2™ cycles mentioned above for which 
the study of the phase ordering in terms of two-state spin 
variables may not be legitimate. 

Consider, for example, a CML with local map 5 M de- 
fined by (||) in a 4-band chaotic regime (p G [pz,P2\) 
which exhibits period-4 NTCB. The "natural" spin vari- 
ables to study phase ordering take four values, indexed 
by the 4-band chaotic cycle. However, these four bands 
can be grouped in two "meta-bands" , since they arise 
from a band splitting bifurcation at p,2 , so that two-state 
spin variables can still be defined. Accordingly, two limit 
coupling strengths can be defined: gl, marking the onset 
of complete phase ordering between the two meta-bands, 
and g 2 for ordering from initial conditions within one of 
the meta-bands. A priori, g 2 ^ gg, and there might exist 
coupling strengths such that, e.g., pinned clusters exist 
within, but not between, the two meta-bands. The "true" 
onset of period-4 NTCB is then given by g > max(g^ , g 2 ). 
Similarly, for p G [poo,Pn], one can define n different p- 
dependent limit coupling strengths gl, g%, ■ ■ ■ , 5™, with 
n — > oo as p — > poo. Using our recent work on renor- 
malisation group arguments for CMLs ||, one can show 
that the threshold values of this infinite hierarchy are re- 
lated to each other. Here, we only describe briefly these 
results, while a detailed derivation can be found in [||. 
The RG structure of single map (||) induces the conju- 
gacy between (A™ o S M ) 2 and A 2 g m o S,( M ), where 
transforms each variable X? by S^, A™ is the diffusive 
operator applied m times, and q(p) — p 2 for coupled 
tent maps. This relation can be shown to imply that 
g 2 (p,m) = g\(q(p),2m). Furthermore, using the fact 
that g2(p,m) decreases with m, one can prove that the 
maximum g e for all n, p, and m is g* — g\ (pi, 1). Thus, 
whenever g > g*, complete ordering occurs for all bands. 

The above results are at odds with the behavior of 
usual models studied in phase ordering problems p5[ . 
But in both cases presented here, the exponent ratio 9/cf> 
seems to take the value expected for the TDGLE model. 
This "weak universality" is reminiscent of similar results 
found recently at the Ising-like critical points shown by 
the same models |lq |. We note, moreover, that, when g 
is increased, <fi approaches 1/2 and 9 reaches values close 
to 9ql- We believe that this tendancy is mostly due to 
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the lattice effects becoming less and less important (al- 
though strict pinning does not occur for g > g c ). We have 
shown recently |8| that, in the continuous-space limit of 
CMLs, the weak coupling regime disappears (g e — > 0), 
together with any pinning effects. One can thus wonder 
whether, in this limit, one recovers more "conventional" 
phase ordering dynamics. 

The continuous limit of CMLs such as those defined 
by is reached when applying the coupling step 

of the dynamics more and more times per iteration, 



oo limit of A" 1 o S. 



In 



g converges to a universal Gaussian kernel 



i.e. when taking the to 
this limit. A 

Af = exp(^-V 2 ) with a coupling range A = ^/2<7m||e|| 
where ||e|| is the lattice spacing, which can thus be chosen 
to scale like 1/y/m so as to keep A constant. We investi- 
gated the phase ordering properties of these CMLs with 
the symmetric local map (|[) for increasing values of to. 
At a qualitative level, the scaling behavior of L(t) and 
p(t) is observed at all to values. Quantitatively, expo- 
nents 9 and 4> vary with to at fixed g. Increasing to, <j) 
seems to converge to 1/2, while — > 9q^: for to = 1 to 
3, we find <j) = 0.467, 0.479, 0.505, and 6 = 0.174, 0.184, 
0.196, from single runs on lattices of linear size 4096 sites. 

Our work provides a quantitative method for determin- 
ing the onset of NTCB in chaotic coupled map lattices. 
It also reveals that the phase-ordering properties of mul- 
tiphase, chaotic CMLs are different from those of most 
models studied traditionally. More work is needed, espe- 
cially at the analytical level, to clarify the origin of the 
non-universality observed and put our numerical results 
on firmer ground, since we cannot completely exclude a 
very slow, unobservable, crossover of the scaling behavior 
observed to that of a more traditional model. Different 
approaches can be suggested. 

A continuous variation of the scaling exponent (f> for 
the characteristic length of domains is not usually ob- 
served, but (at least) two exceptions are known. One is 
the case of coarsening from initial conditions with built- 
in long-range correlations Jl7| , but then the persistence 
probability p{t) does not decrease algebraically with time 
ijLsf . Another situation of possible relevance is the case of 
phase-ordering with an order-parameter-dependent mo- 
bility fl9|| , for which, unfortunately, the behavior of the 
persistence is not known. At any rate, the recovery of 
the "normal" scaling properties of the TDGLE in the 
space-continuous limit suggests that lattice effects are 
ultimately responsible for the non-trivial scaling proper- 
ties recorded in discrete systems. This calls for a detailed 
study of interface dynamics in order to assess the effective 



role of discretization and anisotropy. 

Finally, we believe our results are general and that 
similar behavior should be found in experiments on 
phase-ordering of pattern- forming systems, such as, e.g., 
electro-hydrodynamical convection in liquid crystals, or 
Rayleigh-Benard convection [^o| . 

We thank Ivan Dornic for many fruitful discussions and 
his keen interest in our work. 
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